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In this work we propose a subtle change in Axehod's model for the dissemination of culture. The 
mechanism consists of excluding non-interacting neighbours from the set of neighbours out of which 
an agent is drawn for potential cultural interactions. Although the alteration proposed does not 
alter topologically the configuration space, it yields significant qualitative changes, specifically the 
emergence of surface tension, driving the system in some cases to metastable states. The transient 
behaviour is considerably richer, and cultural regions have become stable leading to the formation 
of different spatio-temporal structures. A new metastable "glassy" phase emerges between the 
globalised phase and the polarised, multicultural phase. 



I. INTRODUCTION 

The interest in phenomena involving information flow 
on networks has been increasing over the last decades. A 
variety of mathematical techniques have been combined 
to model brain activity, biochemical networks, excitable 
media, ecosystems and the like [1]. Naturally amongst 
them are the social systems, including human societies 
and culture. In an attempt to escape from reduction- 
ism, and inspired by the successes of statistical mechan- 
ics, many different approaches are being proposed to un- 
derstand the collective behaviour that emerges from the 
interaction of many subparts that constitute the com- 
plex system under consideration. There is already a con- 
siderable amount of literature on social modelling, from 
applause synchronisation to disease spreading, opinion 
formation and urban spatial segregation. For a general 
review, refer to [2]. In this context, Axelrod proposed in 
1997 a discrete vector model for cultural dissemina- 
tion that has been attracting scientists for its very simple 
description but rich behaviour. 

Axelrod's model fits in a broad class of models based 
on agents, where each agent has an internal state and 
a set of rules determine the mechanisms of interaction 
with other agents. Although very simple in Axelrod's 
case, this framework is very powerful and its possibilities 
of representing multilevel structures and capturing the 
essence of quite different systems are vast. In a certain 
way, Axelrod's model extends some of the basic mod- 
els proposed in the previous decades, such as the voter 
model, to a discrete vector representation of the internal 
state of each site (or node) of the network. Within this 
description, agents can be compared by overlapping their 
vectors. Axelrod's basic premise comes from the fact 
that more similar individuals tend to exchange culturally 
more often, usually referred to as homophily. Though one 
agent gets more similar to its neighbour after the interac- 
tion (local convergence), the global state of the cultural 
network may in some cases tend to a polarised configu- 



ration, where groups of different cultures form. 



Probably one of the most important characteristics of 
Axelrod's model is its out-of-equilibrium phase transi- 
tion, capturing two different regimes: one leading to a 
globalised state and the other to a polarised one, depend- 
ing on the model's parameters. Conversely, not much at- 
tention has been given to its transient behaviour and how 
the spatial structures develop over time. In fact, although 
homophily has been the source of inspiration for Axel- 
rod's mechanics, during time evolution the information 
flows through space - cultures hop through agents quite 
freely - and cultural borders in general tend to break, as 
if people from the same culture did not tend to stick to- 
gether. To account for a stronger homophily we should 
think of a surface tension effect, for which we propose 
a mechanism, that is not directly related to a memory 
or majority rule, but on the contrary, emerges naturally 
from a slight and justifiable change in the updating rules 
of the model. 



In this work we will outline an analysis of the effect 
of introducing a sort of agent optimisation, in the sense 
that agents only use their time to interact with whom 
the interaction is possible. In the next section, we will 
shortly describe Axelrod's original model and some im- 
portant quantities to characterise the macroscopic state 
of the system. In the following section we introduce our 
slight variation in the original dynamics and the moti- 
vation for doing it. The subsequent section contains a 
description of the difference attained both in transient 
and in stationary behaviour as well as some results from 
numeric simulations. Finally in the last section we will 
discuss the resulting differences and present some of the 
new issues that have arisen from the introduction of sur- 
face tension. 
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II. AXELROD'S MODEL 

In the original paper [3 , Axelrod defined culture - 
in a simple fashion - as the set of people's character- 
istics subject to other individuals' influence. In this 
model, every agent i is represented by a discrete vec- 
tor ai = {aj^af^ . . . , a/, . . . , af ) of F components (or 
cultural features), each of which can express one of q 
possible cultural traits (a/ G {0,1,..., (7 — 1}). Each 
feature denotes one cultural aspect, such as musical pref- 
erences, sports, and so forth, and the cultural traits can 
be thought of as different expressions of these features. 
The fact that the number of possible traits per feature 
is the same for all features is just a consideration of the 
model to make it simpler. 

The dynamic of the model is driven by the empirical 
observation that similar individuals are more likely to in- 
teract than dissimilar ones (homophily). The agents are 
situated in the nodes of a network. In this work, we have 
studied the case of a square lattice with periodic bound- 
ary conditions, but other topologies have been used [4], 
[5]. This is a discrete stochastic and asynchronous model. 
In every time step one agent i is randomly selected along 
with one of its neighbours j. The set of neighbours of a 
node i (its neighbourood) is denoted as Ui. An interac- 
tion occurs with probability ooij equal to their similarity 

- the similarity between two vectors ai and aj is calcu- 
lated as their overlap normalised by the number of fea- 
tures Uij = Xl/=i ^cr^ ' where S^f ^/ is the Kronecker 
delta - and it consists of randomly selecting one of the 
features (/) whose trait is different between the agents 
and making the agent j adopt z's trait for that feature 
after the interaction, that is, = a/. It is worth not- 
ing that, in Axelrod's model, it does not matter whether 
i transfers one trait to j or the other way around, it is 
symmetryc with respect to that. 

Thus, two agents i and j are able to interact if their 
similarity lies in the range < cuij < 1, and we will call 
this link active. If ujij = 0, their similarity is null and so 
is their probability of interacting, according to the defi- 
nition. If uJij = 1, there is no more cultural difference be- 
tween them, and no more possible information exchange. 
In both cases, the link between these two agents is said 
to be inactive. One site (agent) is said to be active if it 
has at least one neighbour with whom the link is active. 

Under these rules the system always orbits to an ab- 
sorbing state (without any more possible interactions) 
which can - depending on the two parameters F and q 

- be a globalised or a polarised absorbing state. The 
globalised (monocultural, ordered) state is characterised 
by the majority of agents sharing the same culture (a 
giant component of one culture), whereas the polarised 
(multicultural, disordered) state is characterised by many 
local cultural clusters but no global dominant culture. 
On most topologies, a first order (except for the special 
case F = 2) out-of-equilibrium phase transition separates 
these regimes: there is a critical line in the parameter 



space qc{F) separating these two different phases. Sev- 
eral order parameters have been proposed to characterise 
the transition, but the most widely used is the average 
fraction of the dominant culture's size, i^?^pi. As the 
lattice size increases, the transition gets sharper across 
the critical value (figure [T]), typical of a first-order transi- 
tion. It is worth mentioning that the (out-of-equilibrium) 
transition occurs going from one system to another: q is 
part of the definition of the system, and not a collective 
emergent property such as the temperature. The tran- 
sition happens throughout the absorbing states of these 
different systems as the parameters vary. 
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FIG. 1. ^"^^^^^ vs. q in the absorbing states for several lattice 
sizes. The curves represent L = 5 (black), L = 10 (red), 
L = 20 (green) and L = 40 (blue), averaged over 50, 500, 
1000 and 2000 independent realisations, respectively. F was 
set to 5. 



III. THE MODEL WITH SURFACE TENSION 

The introduction of surface tension in models on net- 
works is sometimes achieved by the use of a majority 
rule [6 or a local memory [7 . In the former case, the 
configuration space topology is altered - some transitions 
otherwise possible are no longer possible with the major- 
ity rule. In this work we propose a subtle modification of 
the model's dynamics that leads to the emergence of sur- 
face tension without changing the topology of the config- 
uration space. All the possible transitions are preserved, 
only the probabilities associated with these transitions 
change. 

The main idea behind this modification is a sort of op- 
timisation from the agent's point of view. Suppose - in 
the Axelrod's case - an agent i was selected and it has 4 
neighbours, one of which (j) has a state such that uJij = 
or UJij = 1. Even though there is no possible interaction, 
there is a finite probability (|) that this neighbour is se- 
lected and the iteration will have to be discarded. People 
usually do not try to interact culturally with whom there 
is no possible interaction, and this is the only alteration 
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introduced in the original model. We can now formalise 
the rules for the model. 

The agents are placed on a network with a specific 
topology (in our case, a 2D regular lattice with periodic 
boundary conditions). The concept of active sites (or 
agents) and active links according to the definition in the 
previous section is still used here. The initial conditions 
are randomly selected from a uniform distribution and 
the dynamics consist of iterating these two steps: 

Step 1: Randomly choose one site i. Then choose one 
of its neighbours j such that < Uij < 1, in case 
there is (otherwise go back to step 1). 

Step 2: With a probability equal to their similarity, 
there is an interaction. The interaction consists 
of randomly selecting an / such that a{ ^ <jj and 

making site j adopt a/. 

Note that the only difference compared to Axelrod's 
original model is the constraint that the randomly se- 
lected neighbour has < uJij < 1. This choice makes 
the model asymmetric with respect to the trait transfer 
direction. If we invert and make site i adopt <jj in step 2, 
the opposite effect arises: the total interface length tends 
to increase and a physically implausible cultural mixing 
occurs. Suppose now a site i is selected and only one of 
the four links connecting it to its neighbours is inactive. 
In this case, there are only three possible neighbours to 
be randomly drawn, with probability | for each. After 
the choice, there is still a probability that no interac- 
tion occurs, but all the selected pairs may potentially 
exchange cultural information. This single modification 
is responsible for the emergence of surface tension (see 
figure |2| . 

A colour scheme was devised for the simulations in 
order to observe the spatiotemporal patterns. Different 
colours were assigned to different cultural states. Being 
the RGB colour space three-dimensional, it is impossi- 
ble to assign similar colour to similar cultures based on 
distances. Thus, the colour to culture mapping does not 
indicate that similar colours represent similar cultures. 

IV. NUMERICAL RESULTS 

While seemingly subtle, there are some radical changes 
in the overall properties of the model. In Axelrod's 
model, the absence of surface tension leads to instability 
of cultural regions, and bubbles tend to break [8 (figure 
[3]). In our model, cultures have a tendency to organise 
in space, forming well-defined cultural regions through- 
out the trajectory in time. This can be interpreted as 
cultural cohesion, or a tendency of people from the same 
culture to stick together: the principle of homophily in a 
stronger sense. As a consequence, instead of local cul- 
tural interfaces there are spatially extended interfaces 
between cultures, macroscopic membranes, which usu- 
ally do not break. These membranes have very peculiar 



FIG. 2. Detail of a square border between two cultures, one 
represented in white and the other in orange. The agents lying 
on the straight part of the interface have the same number 
of active neighbours on both sides (green), thus being the 
probability of trait transfer the same in both directions. In the 
corner, however, that symmetry no longer holds (blue). There 
are two agents which can be drawn to transfer their traits to 
the corner agent, making it twice more probable for the corner 
to be invaded, producing two new corners and continuously 
"rounding" the sharp corners. This illustrates the curvature- 
driven surface tension of the model, not present in Axelrod's 
original model. 

interacting mechanics, and a more careful investigation is 
yet to be made. The isolation of these membranes is pos- 
sible exporting these rules to the voter model, according 
to [9 . Another remarkable consequence of the forma- 
tion of domains is that the cultural exchange happens 
in the borders between regions, and the borders tend to 
concentrate cultural diversity (see figure [4|. 




FIG. 3. Same initial conditions (square), the first row corre- 
sponding to the case with surface tension and the second to 
the Axelrod's original model. 

The transient regime with Axelrod's original dynam- 
ics has in general no tendency for cultural spatial struc- 
tures to form in time. In contrast, after introducing sur- 
face tension it is possible to observe some different spa- 
tiotemporal patterns. For small enough values of q - the 
region of the parameter space where the globalised ab- 
sorbing state is reached - the cultural exchange rapidly 
generates a homogeneous nucleation, covering the whole 
lattice with local polygonal cells whose borders (or mem- 
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FIG. 4. Configuration during transient for (left) F = 11 and 
q — 2^ showing the homogeneous nucleation of bubbles with 
non-rigid membranes and multiculturality on the borders be- 
tween different regions and (right) F — 2b and g = 10, show- 
ing the non-homogeneous nucleation of bubbles, which will 
grow and reach a metastable state. 

branes) are flexible enough, allowing bigger ones to coa- 
lesce progressively, becoming even bigger at the expense 
of anihilating the small ones (figure |4|. For large enough 
values of there are some nuclei with possible culture 
exchange immerse in a vast disordered ocean of cultural 
disagreement. In this region, the transient is very similar 
to Axelrod's, despite the formation of membranes within 
the islands of cultural exchange. The absorbing state 
reached is also similar to Axelrod's in this case. 

Useful to describe the macroscopic state are the densi- 
ties of links, divided in three categories: (inactive) links 
between nodes with no feature in common, (active) links 
between nodes with some features - but not all - in com- 
mon, and (inactive) links between nodes with the same 
cultural state (all features in common). The densities of 
these links will be called po, pa and pp respectively, and 
are given by the following expressions on a regular lattice 
with periodic boundary conditions and N nodes. 

1 ^ 



Pa = ^ - PQ - PF 

One can see how these quantities develop during the 
transient until it reaches a metastable state in the figure 

El 

A function we used to detect metastable configurations 
is the number of surviving traits, proposed in [10 , which 
is a Lyapunov function for both models. It is calculated 
by counting the traits that are still circulating on the 
network, and has maximum value qF and minimum value 
F. Observe that if the numeric value of two different 
cultural features is the same, it is counted twice. The fact 
that they have the same symbol does not mean anything 
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FIG. 5. Graph illustrating the transient behaviour or the 
three categories of bonds: between completely dissimilar sites 
pQ (black), between interacting agents pa (red) and between 
culturally identical agents pp (green), for a lattice size L = 
100, with N = 10000 nodes. In this case, F = 25 and = 
12. The self- averaging makes the curve smoother as lattice 
size grows. A metastable state is reached for this choice of 
parameters (observe that the density of active interfaces does 
not decrease to zero). 

since they stand for different cultural features, thus being 
incomparable. 

/=i t=o 

This function, besides being monotonically decreasing 
in time, may provide insights into the mechanics of the 
models. At each interaction, one agent adopts one of the 
cultural traits of its influencing neighbour and discards 
the previous value for that feature. In this way, there 
are eventual trait extinctions during the system's trajec- 
tory in time. Once there is an extinction of a certain 
trait, it will never get back to the network and many 
configurations previously accessible become inaccessible. 
Kuperman has also studied two different Lyapunov func- 
tions for his model with majority rule [6], and the use 
of these functions to understand stability [11] or detect 
metastable configurations are of great importance. 

What is entirely new is a metastable "glassy" phase for 
intermediate values of q. The transient regime is charac- 
terised by the non-homogeneous nucleation and growth 
of bubbles (figure [4| , and the nucleation usually occurs 
at neighbour sites with orthogonal cultures (in the case 
of non-periodic boundary conditions, the edges and cor- 
ners become nucleation regions). Therefore, the aver- 
age density of nucleation regions also varies with q. Af- 
ter growing, these bubbles eventually touch each other, 
forming well defined domains whose borders are almost 
completely rigid; yet, there might still be some occasional 
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changes. In case of dense nucleation, the mean domain 
size reached is smaher and the domains become more 
susceptible to trait invasion from neighbour ones, pro- 
voking some extra domain concatenation and leading to 
ultimate irregular domain shapes. 

These states are metastable in the sense that they are 
not absorbing states (there are still cultural exchanges 
taking place), but are quasi-stationary, remaining con- 
fined in certain attractor-like regions of the configuration 
space. In [11^ a completely different notion of metasta- 
bility is used, not to be mistaken with the one used here. 
It is important to note that what defines these regions 
as attractors are the transition probabilities, and not the 
topology: if we use these metastable configurations as 
the initial conditions for the Axelrod's case, these do- 
mains break, and the consensus state is reached. Under 
Axelrod's dynamics, there are only fixed point attractors, 
but in our case there are other kinds, characterising a dif- 
ferent attr actor structure, albeit the common topology. 
The absorbing states or metastable states characterising 
the three different phases are represented in figure [6] 




FIG. 6. Absorbing states or metastable states for choices of 
parameters lying in the three different phases. F — 25 and 
L — 500 for all of them. First row: ^ = 2, ^ = 10 and q — 20. 
Second row: q = 50, q = 100 and q — 125. 

The next step would be to characterise both transi- 
tions and find the critical lines in the parameter space. 

This however is not such a simple task. Using the typ- 

(s ) 

ical order parameter ^ "^""^ apparently gives us a good 
initial idea about the critical values but when the lat- 
tice size increases, the middle plateau (corresponding to 
the metastable phase) progressively lowers, tending even- 
tually to zero. We have also tried to use the cluster 
entropy proposed in [12 , which works pretty well for 
Axelrod's original model but not for ours with surface 
tension. One of the difficulties of dealing with this out- 
of-equilibrium transition involving metastable states is 
deciding when to interrupt the simulation to measure 
the order parameter of that (truncated) configuration, 
or even to understand the difference between long tran- 
sients or metastable states. 



V. DISCUSSION 

Axelrod's model is a simplification of a very complex 
problem, and it was not proposed in an attempt to de- 
scribe or explain the intricate process of cultural ex- 
change and formation. That said, the model is inter- 
esting in itself, and that probably is the great motiva- 
tion for studying it. According to our view, one of its 
main limitations is related to the absence of surface ten- 
sion during its time evolution. It is possible to observe 
that cultural traits fiow through the agents in a way that 
the cultural borders tend to break, forming very differ- 
ent spatiotemporal patterns from what we expect due 
to the homophilic nature of social interactions - affine 
agents staying together in space. In Axelrod's model, 
cultural bubbles are however unstable and usually dis- 
solve or break, depending on the surroundings. Besides 
that, both models still lack a mechanism for creation of 
novelty, since the traits only get extinct in the course of 
time. 

We proposed a modification in the rules that results in 
surface tension without changing topologically the config- 
uration space, clearly implying a modification of the tran- 
sient behaviour, making it much richer: cultural borders 
do not break anymore, cultural regions become connected 
in space, spatial patterns emerge such as bubble nucle- 
ation and surface-tension-driven polygonal cells. The 
new model still preserves the globalised and polarised 
phase of Axelrod's model, but a new "glassy" metastable 
phase emerges between them, with cultural regions re- 
sembling political maps whose (active) borders concen- 
trate all the cultural exchange and diversity. Unlike Axel- 
rod's model, this phase has the noteworthy characteristic 
of not reaching a frozen state. This novel outcome ac- 
counts for social interactions in a more realistic way, with 
spontaneous formation of cultural niches and, more im- 
portantly, the boundaries that separate and define them. 

An important issue to be pointed out is the dependence 
of system's evolution on the uniform random distribution 
as the initial condition. The polarised phase actually 
stems from this random initial condition in situations 
where q is large enough such that it is quite improba- 
ble that any two neighbours have something in common. 
Much of our understanding of the model comes from this 
widespread choice, but some carefully chosen initial con- 
ditions lead to significantly different situations that must 
be understood. We have used some of these particular 
configurations as initial conditions to unravel the me- 
chanics of the membranes, and they can be quite useful. 
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